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We investigate some QCD corrections that contribute to the Gross-Llewellyn Smith (GLS) sum 
rule, but have not been included in previous analyses of it. We first review the techniques by 
which the xF% structure function is extracted from combinations of neutrino and antineutrino cross 
sections. Next we investigate corrections to the GLS sum rule, with particular attention to contri- 
butions arising from strange quark distributions and from charge symmetry violating (CSV) parton 
distributions. We find that additional corrections from strange quarks and parton CSV are likely 
to have a small but potentially significant role in decreasing the current discrepancy between the 
experimental and theoretical estimates of the Gross Llewellyn Smith sum rule. 
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I. INTRODUCTION 

The Gross Llewellyn Smith (GLS) sum rule 
(sometimes referred to as the baryon sum rule) 
is obtained from an integral of the xF% structure 
function obtained in charged-current deep inelas- 
tic scattering (DIS) from neutrinos and antineutri- 
nos on nucleon or nuclear targets [U [2] . In recent 
years the most precise neutrino data has been ob- 
tained by the CCFR and NuTeV collaborations, 
from interactions of neutrinos and antineutrinos 
with an iron target From measurements by 
the CCFR group [4] , values have been obtained for 
the GLS sum rule. At the value Q 2 = 3 GeV 2 , the 
CCFR analysis claimed a precision of roughly 3%. 
The CCFR data could also be used to obtain the 
GLS sum rule as a function of Q 2 ; this allows one 
to test contributions from higher-order QCD cor- 
rections [5] and from higher- twist terms [5J. The 
current status of various DIS sum rules has been 
summarized in a review article by Hinchliffe and 
Kwiatkowski [2]. 

In this paper, we point out that some ad- 
ditional QCD effects, particularly contributions 
from strange quarks and parton charge symmetry 
violation (CSV), have not been included to date 
in estimates of the GLS sum rule. One now has 
recent experimental data on strange quark par- 
ton distributions, in particular on the asymmetry 
between strange and antistrange quarks, which is 
relevant for the GLS sum rule. In addition, there 
is much interest in the possibility of CSV in the 
parton distributions [TUTU]. 

We examine the possible contributions of these 



terms to the Gross-Llewellyn Smith sum rule, 
finding that such corrections are likely to pro- 
duce small but potentially significant effects. At 
present, theoretical estimates of the GLS sum 
rule lie one or two standard deviations below the 
data [2J. We point out that contributions from 
strange quarks and partonic CSV are likely to im- 
prove the agreement between theory and experi- 
ment. 

Our paper is organized as follows. In Sect. [IT] wc 
review the form of neutrino cross sections and the 
derivation of the GLS sum rule. The experimen- 
tal results of the CCFR group [3] are summarized 
and compared with the theoretical calculations of 
Hinchliffe and Kwiatkowski [2J. 
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In Sect. 

review how the structure functions, particularly 
xF 3 , are extracted from experimental data. We 
pay special attention to contributions from strange 



quarks and partonic CSV. In Sect. |IV[ we make 
estimates of these contributions to the GLS sum 
rule. 

In Sect. |IV D| we review isoscalar corrections 
to the data, which arise because iron is not an 
isoscalar target. Because of the way in which 
isoscalar corrections have been implemented in 
previous analyses, it is difficult for us to give a 
definitive, quantitative estimate of the contribu- 
tion of strange quarks and partonic charge symme- 
try violation in the GLS sum rule. Nevertheless, 
our analysis clearly establishes that these correc- 
tions can be as large as the quoted errors on the 
GLS sum rule and that they tend to improve the 
agreement between theory and experiment. Our 
results show that these corrections should be in- 



eluded in future analyses. 



II. NEUTRINO CROSS SECTIONS AND 
THE GROSS-LLEWELLYN SMITH SUM 
RULE 

The cross section for charged current (CC) in- 
teractions initiated by neutrinos or antineutrinos 
on nucleons on a proton is shown schematically in 
Fig. [I] It has the form [3 [7] 
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The relativistic invariants in Eq. Q are Q 2 = —q 2 , 
the square of the four momentum transfer for the 
reaction, x and y. For four momentum k (p) for 
the initial state lepton (nucleon), we have the re- 
lations 

x = ^; j/=^|; s=(k+p) 2 . (2) 
Ip ■ q p ■ k 

Eq. {l} applies in the limit Q 2 « M^. We 
have introduced the Fermi coupling constant, G F , 
in terms of the electromagnetic coupling constant 
a, the W boson mass M w , and the weak mixing 
angle 9 W , 



G, 



V2sin 2 9 w Ml 
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to a W + being absorbed on the proton. In Eq. |T]) 
we have written the structure function F\ in terms 
of the structure function F 2 and the longitudinal 



to transverse cross section ratio R v Ll i.e. 
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For the CC reactions initiated by neutrinos, the 
experimental values for R V L are summarized in the 
review article by Conrad, Shaevitz and Bolton [3]. 

The most accurate neutrino and antineutrino 
cross sections are on nuclear targets, particularly 
(in the case of the CCFR and NuTeV measure- 
ments) on iron targets. The structure functions 
are typically described in terms of nuclear parton 
distributions. At sufficiently high values of Q 2 , the 
structure functions per nucleon for a nucleus with 
Z protons and N = A — Z neutrons can be written 
in terms of averages and differences of structure 
functions for neutrinos and antineutrinos. 

Therefore we define 
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FIG. 1: Schematic picture of charged-current ampli- 
tudes in DIS induced by neutrinos or antineutrinos. 
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For neutrino- induced interactions, F 2 
represents the F 2 structure function corresponding 



In terms of parton distribution functions it is 
straightforward to show that 



F \ VA {x) = x(u + (x) + d + (x) + s + (x) + c + (x)) 
- x^(Su + (x) + Sd+{x)) ; 



xF^ (x) = x(u (x) + d (x) + s (x) + c (x)) 
N 

— x — (Su (x) + Sd (x)) ; 

AF 2 WA (x) = 2x[s-{x)-c-{x) + v f (u-(x) 
N 

- d (x)) + —(-5u (x) + 8d~ (x))] \ 
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AxF^ A (x) = 2x[s + {x)-c + {x) + v s {u + {x) 

- d + (x)) + —(-Su + (x)+5d + {x))] . 

(6) 

In Eq. Q we have assumed the impulse approx- 
imation, i.e. that the nuclear structure functions 
are simply given as the sum of free nucleon parton 
distributions. Nuclear modifications of the par- 
ton distributions have been considered by various 
groups [TTHTt)] . and we note particularly the re- 
cent discovery of an isovector EMC effect [17] . In 
Eq. (Jsj) we use the notation 

q ± (x) = q(x)±q(x) . (7) 

We have introduced the neutron asymmetry pa- 
rameter Vf = (N - Z)/A. For the CCFR and 
NuTcV iron targets the average value is Vf = 
0.0567 [HJ Unj. We also include possible parton 
charge symmetry violation (CSV), with the nota- 
tion 
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and an analogous equation for antiquarks. 

The neutrino/antineutrino average structure 
functions contain contributions from light quarks 
including the relatively large light valence quark 
contributions. The structure function differences 
on a nuclear target contain contributions from 
heavy quarks, partonic CSV contributions, and 
light quark contributions proportional to the neu- 
tron asymmetry Vf. In previous work the term 
AxF^ A (x) has been treated as a small perturba- 
tion, while the term AF^ A {x) has been neglected. 
In this paper we will examine the possible effects of 
the term AF^ A (x) on the Gross-Llewellyn Smith 
sum rule. 



A. The Gross-Llewellyn Smith Sum Rule 

The Gross-Llewellyn Smith (GLS) Sum Rule 
[Tj is derived from the first moment of the F3 
structure functions for neutrinos and antineutri- 
nos. The easiest derivation of the Gross-Llewellyn 
Smith sum rule results from summing the xF 3 
structure functions for neutrinos and antineutri- 
nos on a proton. Using Eq. ^ for the case of a 
proton (Z = 1, N = 0) we obtain 
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In Eq. Q we have neglected various QCD cor- 
rections and higher twist contributions. Without 
these corrections, the GLS sum rule is equal to 
three because the first moment of the light quark 
parton distributions gives the total number of va- 
lence quarks in the nucleon. From valence quark 
normalization, the first moments of the strange 
and charm quark asymmetries, s~ (x) and (x) 
respectively, vanish when integrated over all x. 

In order to compare the Gross-Llewellyn Smith 
sum rule with experimental data, we must consider 
a number of corrections. For pedagogical purposes 
we will discuss the GLS sum rule for an isoscalar 
target. At the end of this article we will review the 
corrections that are made for non-isoscalar targets 
such as iron. After applying a series of QCD cor- 
rections one obtains 
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The naive Gross-Llewellyn Smith sum rule is cor- 
rect only in leading twist approximation, and only 
to lowest order in the strong coupling constant 
a s . Our expression for the GLS sum rule thus 
includes a QCD correction (the term in square 
brackets in Eq. (10)), which was derived by Larin 
and Vermaseren [5] using a QCD scale parameter 
A QCD = 213 ± 50 MeV, and the quantity AHT 
represents a higher twist contribution [6 . This 
is summarized in the review article on QCD sum 
rules by Hinchliffe and Kwiatkowski [2]. 

As is the case for the Adler [5D] and Gottfried 
[2"T] sum rules, the Gross-Llewellyn Smith sum rule 
requires that the structure function be divided by 
x in performing the integral. This gives a strong 
weighting to the small- x region, such that as much 
as 90% of the sum rule comes from the region 
x < 0.1. The most precise value has been ob- 
tained by the CCFR collaboration [3] , which mea- 
sured neutrino and antineutrino cross sections on 
an iron target, using the quadrupole triplet beam 
(QTB) at Fermilab. A summary of experimen- 
tal details for precision measurements using high- 
energy neutrino beams is given in the review arti- 
cle by Conrad, Shaevitz and Bolton [3] , and a de- 
tailed description of the experimental details and 
analysis procedure used by the CCFR collabora- 
tion is given in the thesis of Seligman [IS] . 
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Because of the large contribution to the GLS 
sum rule from small x, one measures xF% at vari- 
ous values of x, and evaluates the integral 

W»)= f ldjJ yFY A (y)- (ID 

Jx v 

The Gross-Llewellyn Smith sum rule is then ob- 
tained by taking the limit 

S GLS = lim^o S GLS (x) . (12) 



contributions. We will estimate each correction 
term as a function of x. The lowest x value con- 
tributing to the Gross-Llewellyn Smith sum rule 
as measured by the CCFR group is x min = 0.015. 
We will calculate each contribution to the GLS 
sum rule as a function of x, and estimate the con- 
tribution SS G Ls(x)\ x=Xmin . 
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FIG. 2: Experimental results for Gross-Llewellyn 
Smith sum rule, Eq. [lOj from the CCFR group, Ref. 
[4]. Squares: xFs(x), sum of neutrinos plus antineu- 
trinos, at Q 2 — 3 GeV 2 . Dashed curve: analytic 
fit to XF3. Diamonds: approximation to the inte- 
gral Sgls(x) of Eq. 
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Solid line: fit to the integral 



Fig. [2] shows the CCFR measurements on iron 

and the experimental values of xF 3 (x) (the sum 
of the nuclear xF% structure function for neutri- 
nos plus that for antincutrinos) vs. x. The CCFR 
group measured cross sections at several values of 
x and Q 2 . The squares give the value of xF 3 (x) in- 
terpolated to an average momentum transfer Q 2 = 
3 GeV 2 (this is the mean Q 2 for the lowest x-bin 
in the CCFR experiment, since the lowest x values 
contribute the greatest amount to the GLS sum 
rule). The dashed curve is the best fit to xF^{x) 
of the form Ax b {l — x) c . This form was used to 
extrapolate the first moment to x = 0. The CCFR 
reported value for the sum rule H] at this Q 2 value 
is S G ls = 2.50 ± 0.018 (stat) ± 0.078 {syst). The 
GLS sum rule is therefore known to about 3%. 

The solid curve in Fig. [2] is S GLS (x). In the fol- 
lowing sections we will consider additional QCD 



CHARM 2.56 ± 0.43 

CCFRJ? Z.B3 ± 0.20 

2 70 ± 0.40 
CCFR(NBB) E.7B ± 045 
CCFR(QTB) 2.50 ± 0.0B 



FIG. 3: Gross-Llewellyn Smith sum rule, and errors, 
for a series of experiments, in chronological order from 
top to bottom. 



Fig. [3] shows the evolution over time of the GLS 
sum rule value. The measurements shown are 
from the CDHS [22], CHARM [23], CCFRR [24], 
and WA25 [3S] collaborations. There are also two 
points from the CCFR measurements, the first us- 
ing the Narrow Band Beam (NBB) neutrino data 
[26] 127] and the second using the QTB data [4] 
from the Fermilab Tevatron. 

The points with error bars in Fig. [4] represent 
the experimental results from the CCFR group 
for the Gross-Llewellyn Smith sum rule as a func- 
tion of Q 2 . The curves are theoretical QCD pre- 
dictions by Hinchliffe and Kwiatkowski [2], using 
higher-order QCD corrections from Larin and Ver- 
maseren [S] , and higher-twist corrections of Braun 
and Kolesnichenko [6 . The dashed curves are 
calculations without higher-twist effects, and the 
solid curves include higher- twist. The theoretical 
calculations appear to lie systematically below the 
experimental results by one to two standard devi- 
ations. 

In the remainder of this paper we will review the 
steps that are taken to extract the F 3 structure 
functions from the experimental cross sections. 
We will then review the corrections to Eq. In 
particular, we will focus on the contributions to 
the Gross-Llewellyn Smith sum rule from strange 
quarks and from charge symmetry violating con- 
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FIG. 4: Results for the Gross-Llewellyn Smith sum 
rule, S G ls for various values of Q 2 \ from Ref. 2 . The 
inner error bars are statistical and the outer errors 
combine statistical and systematic errors. The data 
are from the CCFR experiment [4], and the solid and 
dashed curves are theoreteical QCD predictions from 
Hinchliffe and Kwiatkowski [2] ; the solid curve includes 
higher-twist effects while the dashed curve neglects 
them. 



tributions to parton distribution functions. Al- 
though these constitute fairly small corrections, 
nevertheless they may play an important role in 
determining the extent of the agreement, or dis- 
agreement, between theory and experiment for the 
GLS sum rule. 
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In Eq. Q the coefficients f{y,Q 2 ) and g{y,Q 2 ) 
are defined as 
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The cross sections entering Eq. ( 13 1 have been nor 



AF WA 



malized to give the correct total cross sections for 
neutrinos and antineutrinos. The procedure for 
this is described in the review by Conrad et al. [3] . 
For the time being, we will consider the extraction 
of the structure function xF 3 WA for an isoscalar 
target. 

Previous analyses have neglected the term 
(x) in Eq. (13 1. From Eq. ^ we see that 
all contributions to this term should be small. 
We will neglect the c~ term in Eq. (|9| because, 
even though a mechanism has been identified 
which could produce such an asymmetry |28j . the 
charm contribution is certainly suppressed sub- 
stantially with respect to that associated with 
strange quarks and the charmed contribution is 
also kinematically suppressed for the experimcn 
tal conditions 



c.f. Eq. (37) in Sect 



for an isoscalar target at sufficiently 
expect 
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righ Q 2 we 



III. EXTRACTING STRUCTURE 
FUNCTIONS FROM NEUTRINO AND 
ANTINEUTRINO CROSS SECTIONS 
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The most precise neutrino and antineutrino DIS 
cross sections are those of the CCFR and NuTeV 
groups, both taken on iron targets. From Eq. (JlJ 
by taking sums and differences of differential cross 
sections for charged-current DIS from neutrino 
and antineutrino beams, we can isolate different 
combinations of structure functions. The CCFR 
and NuTeV experiments bin the data in x and 
Q 2 . Defining the quantity c = n/G^ME^, it is 
straightforward to show that 
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From Eq. ( 15 ) we see that for an isoscalar target 



the term AF 2 A {x) will be non-zero only if one has 
a strange quark momentum asymmetry s~(x) ^ 0, 
and/or non-zero valence quark CSV contributions. 
In Sect. lIVDl we will discuss additional contribu- 
tions to AF 2 WA {x) for a non-isoscalar target. 
If one neglects the term AF 2 WA (x) in Eq. (Tt3} 



(as was the case in the analysis of the CCFR data), 

then the structure function xF 3 (x) will just be 
proportional to the difference between the neu- 
trino and antineutrino charged-current DIS cross 
sections. For a given x bin, the structure function 
xF 3 is then given by averaging the structure func- 
tion differences over the Q 2 bin appropriate to the 
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given x bin. Thus we obtain 
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In Eq. (16), (Q ) denotes the average over the Q 



bin appropriate to a given x bin. In Eq. (16 1, we 

— wa — 9 

have neglected the slow variation of F 3 with Q . 
However, if the quantity AF^ A (x) is non-zero, 



then Eq. (16) will not give the structure function 
xF 3 (x), but rather a linear combination of xF 3 



and AF2. Comparing this with Eq. (13), we note 



that the y dependence of the coefficients / and 
g of the two terms is quite different, as shown in 

I I — WA 

Eq. (14). In particular, the coefficient of xF 3 
vanishes at y = 0, while the coefficient of AF^ A 
is finite. Inserting Eq. ( 13 ) into Eq. ( 16 ) and av- 



eraging over the Q range for each x bin gives 
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The quantity B(x) in Eq. (17) is the relative 
weighting between the AF 2 and xF 3 terms. B(x) 
will depend upon the a; bin and the Q 2 values that 
are averaged over for each x bin. In Fig. [5] we 
plot the ratio f(y,Q 2 )/g(y,Q 2 ) in Eq. |l3| vs. y 
(this quantity is identical to the ratio fi(y)/f2(y) 
from Eq. ([I])). This ratio is always greater than 
one, and becomes quite large at small y values. 
The quantity B(x) is normalized to equal one if 
one integrates over all y; however B{x) could be 
greater than one particularly if the average over 
Q 2 is weighted towards small y values. 

Eq. (17) shows that for an isoscalar target, the 



common process of taking the difference of neu- 
trino and antineutrino charged-current cross sec- 
tions, and averaging over the Q 2 bin for a given 

— WA 

x bin produces a linear combination of xF 3 
and AF^ A with a relative weighting B{x). From 
Eq. (|6| we see that the partonic content of the 
quantity identified as xF 3 will be 
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xF 3 "(x)" = xF 3 "(x) + B{x)AFY a {x) 
x[u~{x) + d~(x) + s~(x)(l + 2B(x)) 

Sd-(x)(B(x)- 1 -)~6u-(x)(B(x)+ 1 -)]. 

(18) 




FIG. 5: The quantity f(y,Q 2 )/g{y,Q 2 ) from Eq. |l3| 
plotted versus y, under the assumption that W L = 0; 
it is identical to /i(j/)//2(y) defined in Eq. |l]). 



— WA 

We put quotes around the quantity xF 3 , since 

— WA 

it represents the linear combination of xF 3 and 
AF% VA obtained for that x bin. As we have 
mentioned, the y dependence of the coefficients 
/(2/iQ 2 ) an d 9(y,Q 2 ) in Eq- ( 13 1 is quite differ- 
ent, particularly in the forward direction. When 
the longitudinal to transverse ratio R V L = 0, the 
ratio f(y,Q 2 )/g(y,Q 2 ) is given in Fig. [5] If we 
assume that for each x bin, the data is averaged 
over all y, then one obtains B(x) — 1, and Eq. ( 18 ) 
becomes 
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3 ( X ) \B(x)=l 
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(19) 



These additional terms should be included in the 
analysis of the experimental data. In the absence 
of such an analysis we will provide estimates of the 
sign and magnitude of these corrections and their 
effect on the Gross-Llewellyn Smith sum rule. 



IV. ADDITIONAL CORRECTIONS TO 
THE GLS SUM RULE 

In addition to contributions from the light va- 
lence quarks, higher order QCD terms and higher 
twist contributions, Eq. ( 10 ) contains additional 
QCD corrections. These corrections are included 
in the experimental determination of the GLS sum 
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rule, but these additional terms have not been in- 
cluded in theoretical calculations. There are two 
types of corrections. The first involves additional 

contributions to the desired term xF 3 in Eq. (161) . 
The second contribution appears by virtue of the 
term AF^ A , which has not been separated from 
the desired term. In both cases we include con- 
tributions from strange quarks and partonic CSV 
corrections. Although the first moment of each 
of these contributions vanishes when integrated 
over all x, there is a residual contribution still 
present at x m i n — 0.015, the lowest data point 
in the CCFR experiment. In order to reconcile 
theoretical and experimental determination of the 
GLS sum rule, we choose to subtract these ad- 
ditional contributions from the experimental data 
when comparing with theory. We define the cor- 
rections to the GLS sum rule as 



the C-odd strange quark and valence CSV effects 
contribute to SS GLS {x) at any finite value of x. So 
even if the quantity B(y) were a constant, the CSV 
and strange quark effects would be finite for any 
non-zero value of x, vanishing only at x = 0. 

From our current understanding of the parton 
distributions, for sufficiently large x we expect ev- 
ery term in the quantity AF^ A (x) in Eq. (15 1 to 
have the same sign. As we shall see, all of the lat- 
est analyses of strange quark distributions [29TI33] 
find that the quantity s~ (x) is positive for suffi- 
ciently large x. Analyses of parton valence charge 
symmetry violating effects for parton distributions 
[7] obtain a quantity 6d~ (x) — 8vT (x) > for 
x > 0.1. Consequently, both of the terms in 8S GLS 
in Eq. (21 1 will contribute with the same sign. In 



the following sections we will estimate the magni- 
tude of each of these contributions. 
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For an isoscalar target these corrections have the 
form 

SS GLS (x) = f(a s )[5S GLS (x) + SS^(x)} , 
SS s GLS (x) = f s-(y)(2B(y) + l)dy , 



(B(y)--)Sd-(y) 



(B(y) + -)6u-(y) 



dy 



(21) 



The term f(a s ) in Eq. (21 ) is the QCD correction 
that has been calculated by Larin and Vermaseren 
5 ; it is the term in square brackets in Eq. (10). 
Eq. (21) contains two terms. The first is the 



contribution from the strange quark asymmetry. 
The second is the contribution from charge sym- 
metry violating valence PDFs. An additional ef- 
fect will result from the nuclear modification of the 
parton distributions. Implicitly, all of the parton 
distribution functions in Eq. (plj) denote parton 
distributions in iron. In Sect. HVDI wc will discuss 
nuclear modifications of the PDFs. Note that the 
terms containing the quantity B(y) result from the 
contamination of the xF 3 structure function from 
the AF2 term. 



Corrections to the GLS Sum Rule at 

Q 2 = 8 GeV 2 



For convenience we calculate the corrections as- 
sociated with a strange quark asymmetry and par- 
ton charge symmetry violation at a single value of 
Q 2 , and we choose Q 2 — 8 GeV 2 . From Fig. [4] we 
see that this is a value of Q 2 for which the exper- 
imental GLS sum rule has been measured. It is 
also a relatively convenient value of Q 2 for which 
to estimate corrections from both strange quarks 
and parton CSV. At Q 2 = 8 GeV 2 , the experimen- 
tal value of the Gross-Llewellyn Smith sum rule is 
S% X ls = 2.76 ± 0.14, and the theoretical value of 
Hinchliffc and Kwiatkowski [5] including higher- 
twist corrections is S GLS = 2.62. So the theoreti- 
cal value is just over one standard deviation below 
the experimental result. 

In the next sections we will provide estimates for 
the strange quark and partonic CSV contributions 
to the GLS sum rule. After evaluating them we 
will determine the correction that they make to 
the experimental value and error for the Gross- 
Llewellyn Smith sum rule for this value of Q 2 . 

The QCD correction that appears in Eq. (21) 



was evaluated by Larin and Vermaseren [5] ; it has 
the form 



If the quantity B(y) in Eq. (plj) was a constant, f(a s ) = 1 - — - (—) ri-f—") r 2 , (22) 

mi 1-tntti tho stramrp anH (~NV tprms wmilH iri-up 7T V 7T / \ 7T / 



then both the strange and CSV terms would give 
zero in the limit x — > 0. This is because va- 
lence quark normalization requires that the va- 
lence strange quark and valence CSV PDFs have 
zero first moment. However, we have no reason to 
believe that B(y) will be constant. Note also that 



where for nf — 3 active flavors one has r\ = 
3.5833 and r 2 = 18.9757 0. For the strong 
coupling a s we use the value chosen by Hinch- 
liffe and Kwiatkowski [5] in their theoretical cal- 
culations. In the MS factorization scheme they 
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(51 

chose the scale parameter = 233 MeV; here 

the superscript denotes nj. This corresponds to 
A^^ = 366 MeV and a strong coupling a s (Q 2 — 
3 GeV 2 ) = 0.336. Using this value we then obtain 
f(a s )(Q 2 = 8 GeV 2 ) = 0.885. 



B. Contributions from Strange Quarks 

The strange quark parton distributions are best 
obtained from an analysis of opposite-sign dimuon 
production in reactions induced by neutrinos and 
antincutrinos. In such reactions, dimuon produc- 
tion from a v (y) beam is sensitive to the s (s) 
distribution, so that in principle a comparison of 
these cross sections could enable one to determine 
differences between the s and 5 PDFs. There 
are recent measurements of these reactions by the 
CCFR and NuTeV HQ [35] collaborations. In 
the CCFR experiment the v and V beams were 
not separated and the type of reaction was in- 
ferred from the charge of the faster muon, while 
the NuTeV experiment used separated v and V 
beams. The correction to the GLS sum rule is ob- 
tained from 




0.05 0.1 0.15 0.2 0.25 0.3 0J5 0.4 0.45 



FIG. 6: [color online] The quantity xs (x) = x[s(x) — 
s(x)], vs. x, as extracted by the NuTeV Collaboration, 
Ref. [29]. Values are obtained for Q 2 = 16 GeV 2 . The 
outer error band is the combined error, while the inner 
band is without the uncertainty in the semileptonic 
branching ratio B c . 



(y)(2B(y) + l)dy (23) 



Now, the first moment of is zero, from 
valence quark normalization (there are no net 
"strange valence" quarks in the nucleon). How- 
ever, recent phenomcnological analyses of strange 
quark distributions all obtain qualitatively similar 
results. All of them find the most probable value 
is a positive strange quark momentum asymmetry, 
(xs~(x)) > 0. Also, the best fit to the quantity 
s~ (x) changes sign at an extremely small value of 
x and is large and positive down to rather small x 
values. 

For example, the analysis by Mason et al. [29] 
obtains a best value for the the integral of xs" (x) 

S- = 0.00196 ± 0.00046(stai) ± 0.00045(syst) 
±° Zt*(external). (24) 

The quantity Q refers to the second moment of a 
parton distribution q(x), i.e. 



Q 



x q(x) dx 



(25) 



The Mason result is obtained for a value Q 2 = 
16 GeV 2 . In Eq. (24 1, the term "external" refers 



to the contribution arising from uncertainties on 
external measurements. 

Fig. [6] plots the quantity xs~(x) vs. x from the 
latest NuTeV analysis. The strange quark mo- 
mentum asymmetry, S~ , is quite sensitive to two 
quantities. The first is the semileptonic branching 
ratio B c ; the outer band in Fig. [6] shows the result 
for 5 _ with the B c uncertainty, and the inner band 
is the result without the B c uncertainty. The sec- 
ond is the point at which the distribution xs" (x) 
crosses zero (it must cross zero at least once to give 
zero first moment for s~(x)). The current best fit 
crosses zero at a very small value x ~ 0.004. This 
means that the quantity s~ (x) would have a large 
negative spike at extremely low x (in fact, an x 
value smaller than the lowest x point measured in 
the experiment). 

There are several other recent phenomenolog- 
ical estimates of the strange quark asymmetry. 
All are very sensitive to the CCFR and NuTeV 
dimuon production data. The CTEQ group [30] 
obtains S~ = O.OOlSI^o^ at their starting scale 
Ql = 1.69 GeV 2 . Their best fits to xs~(x) found 
crossovers in the vicinity xq ~ 0.01 — 0.02. The 
NNPDF Collaboration [3T] used only the NuTeV 
data and not the CCFR results. They report a 
value S- = 0.0005 ± 0.0086 at Q 2 = 20 GeV 2 ; 
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S- = {xs-} xo g 2 (GeV 2 ) 

Mason et al. [29] 0.00196 ± 0.00143 0.004 16 

CTEQ [30] 0.0018l£;g^ 0.01 - 0.02 1.69 

NNPDF [31] 0.0005 ± 0.0086 0.13 20 

MSTW [32] 0.0016115 ^9 0.016 1.0 

Alekhin et al. [33] 0.0013 ± 0.0009 ± 0.0002 < 0.02 8 



TABLE I: A summary of recent phenomenological estimates of the strangeness asymmetry (S ), the crossover 
point xo and the Q 2 value at which the asymmetry is obtained. 



the exceptionally large error in the NNPDF value 
(a factor of 5 to 6 larger than the errors from the 
other analyses) results in part from their use of 
a neural network procedure, which does not build 
in widely accepted constraints on the shape of sea 
quark distributions. 

Two other groups supplemented the NuTeV and 
CCFR dimuon data with charm production data 
from CHORUS [36l [37] , which helps to constrain 
the semileptonic branching ratio. The MSTW 



group [32] obtains S~ = 0.0016: 



0.0009 



at Q 2 



10 



GeV 2 ; their next to leading order (NLO) fit had a 
crossover xo = 0.016 at the starting scale Ql = 1 
GeV 2 . Alekhin, Kulagin and Petti [33J obtained 
S~ = 0.0013 ± 0.0009 ± 0.0002 at Q 2 = 8 GeV 2 , 
and their best fit to xs~(x) had a crossover xq < 
0.02. For these various phenomenological fits we 
summarize the values of S~ , the crossover point 
Xq and the value of Q 2 at which the asymmetry is 
calculated in Table I. 



From Eq. ( 23 ) the contribution to the GLS sum 



rule from strange quarks, for a given value of x, 
will be given by the integral of s~(y) weighted by 
the quantity 2B(y) + 1. If we approximate B{y) 
as a constant we just need the integral of s~{y) 
from x to 1. We choose the smallest value x m i n = 
0.015 measured by CCFR. To estimate the integral 
we made an analytic fit to the strange asymmetry 
measured by Mason et al., in the region x > 0.004 
[29] . Our fit had the form s~(x) — ax b e~ cx (x — 
0.004). With this fit we obtain the result 



s~{y)dyK 0.026. 



(26) 



x=0.015 



Within their error bars, all of the phenomenologi- 
cal fits now obtain a positive value for the strange 
quark momentum asymmetry. Because of their 
unusually large error on the strange quark asym- 
metry and its unphysically large value in the va- 
lence region, we do not use the NNPDF result. 
All of the other fits to the strange quark PDFs 



produce a quantity xs~(x) which changes sign at 
an extremely small value of x. Thus all of these 
strange quark PDFs will give a reasonably large 
contribution to the GLS sum rule at a value of x 
corresponding to the lowest x value measured in 
the CCFR experiment. We assign an error of 75% 
which is roughly the average of the four determi- 
nations (excluding NNPDF) summarized in Table 
I. Thus we choose 



s~(x)dx = 0.031 ±0.023. 



(27) 



x=0.015 



In Eq. ( 27 1 we have increased the integral of this 
distribution by 20%; this represents the approx- 
imate increase in this moment in evolving from 
Q 2 = 16 GeV 2 to the value Q 2 = 8 GeV 2 ap- 
propriate for our evaluation of the GLS sum rule. 
This increase is comparable to results obtained by 
the NNPDF group [31], who performed DGLAP 
evolution on the second moment of strange quark 
distributions to extrapolate in Q 2 . 

From Eq. ( 23 ) , the strange quark asymmetry 



contribution to the GLS sum rule will depend on 
estimates of B(x). We make two simple guesses for 
this quantity. First, we take B{x) — 1; this is the 
result if the integral of Eq. ( 17) was taken over all 
y. Next, we estimate B(x) = 2; this represents an 
upper limit (and possibly an overestimate) of this 
quantity. Under these approximations we obtain 
strange corrections to the GLS sum rule 

SS S GLS = 0.093 ±0.063 , for B{x) = 1 ; 
5S S GLS = 0.156 ±0.115 , for B{x) = 2 . (28) 

The large result for the strange quark contribution 
results from the fact that the strange quark mo- 
mentum asymmetry changes sign at an extremely 
small value of x. Although we have used the re- 
sult of Mason et al., this property is shared by all 
phenomenological strange quark analyses in Table 
I except for NNPDF. 
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The contribution from the strange momentum 
asymmetry to the GLS sum rule is strongly de- 
pendent on the crossover point xq at which s~ (x) 
crosses zero. If the crossover point for s~(x) oc- 
curred at a value xq > 0.1 then strange quarks 
would make an extremely small contribution to the 
GLS sum rule. It is difficult to imagine a physi- 
cal mechanism that would cause s~ (x) to change 
sign at such small crossover points xq as have 
been found in these phenomenological analyses 
[2"9l |3"UI I3"2l I3"3"] . Indeed, model calculations almost 
invariably yield a zero at x ~ 0.1 or higher [5BTI4T] . 
The NuTeV group found that with a moderate in- 
crease in x 2 one could obtain considerably larger 
values of Xq and corresponding large decreases in 
the second moment S~ 129 . 



C. Contributions from Charge Symmetry 
Violating PDFs 

We also can estimate the contribution from par- 
ton CSV. For an isoscalar target this is given by 



[{B{y)-\)5d-{y) 



(B(y) + -)6u-(y)]dy (29) 



For this we need the valence CSV parton distribu- 
tion functions [7J. We adopt the functional form 
used by the MRST group [3], 

5vT(x) = -Sd~(x) = Kx-°- r °{l - x) 4 (x- .0909) 

(30) 

The best fit value of MRST was k = -0.2. This 
produced contributions very much like the quark 
model valence CSV calculations of Rodionov et 
al. 02] evaluated at Q 2 = 10 GeV 2 . Here we 
choose k — —0.3, which approximates quite well 
the quark model valence CSV from Rodionov, plus 
the valence CSV arising from "QED splitting" 
[43l |44]. If we assume that Sd~(x) = —5u~(x), 



has the form 



then from Eq. (10 1 the term Sg™(x) in Eq. (29 I 



B(y)Su-(y)dy . (31) 



We insert the analytic form of Eq. (30) into 
Eq. (31) and evaluate at x m i n — 0.015, the mini- 



mum x value for the CCFR measurements. As for 
the strange quark contribution we evaluate this us- 
ing two different values for the weighting function, 



B(y) — 1 and B(y) — 2. We assign a 100% error 
to the CSV contribution. Thus we obtain 

SSats = 0.013 ± 0.013 ,B(y) =1 

= 0.026 ±0.026 ,B(y) = 2 (32) 

The partonic charge symmetry violating contri- 
butions correspond to a value Q 2 ~ 10 GeV 2 . This 
is sufficiently close to the value Q 2 = 8 GeV 2 that 
we do not modify this further. 



D. Non-Isoscalar and Nuclear Corrections 

Until now our equations have assumed an 
isoscalar target. We must include corrections to 
account for the excess neutrons in iron. From 
Eq. ^ there are a number of small corrections 
in the case that N ^ Z. For our purposes the 
most important will be an additional contribution 
to the quantity A.F^ VA {x) of the form 



5(AF™ A (x)) 



Vfx[u (x) — d (x)] 



(33) 



When multiplied by the quantity B(x), divided 
by x and integrated over all x, this leads to the 
contribution 

SS%. 9 {x) = v f [ B(y)[u-(y)-d-(y)}dy. (34) 



For convenience we take the lower limit of this in- 
tegral to be x = 0; in the limit where B(x) is ap- 
proximated as a constant, B{x) ~ B, the remain- 
ing integral is just one, so in this approximation 
we obtain 



SScis ~ VfB , 
= 0.058, 
= 0.115, 



for B(x) = 1; 
for B(x) = 2. 



(35) 



Note that all of the contributions to the GLS 
sum rule (strange quark asymmetry, quark model 
valence CSV and QED splitting CSV, and non- 
isoscalar effects) are of the same sign. Thus, their 
contributions will add coherently in modifying the 
GLS sum rule. 

The contributions listed here for strange quarks, 
CSV parton distributions and N ^ Z effects 
were included in the experimental results but not 
in the theoretical calculations of Hinchliffc and 
Kwiatkowski [2]. Consequently we choose to sub- 
tract these contributions from the experimental 
results, in order to compare with theory. At 
Q 2 = 8 GeV 2 the quoted experimental result for 
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the Gross-Llewellyn Smith sum rule was 2.76 ± 
0.14. We take the strange quark contribution from 
Eq. ([28]), the CSV contribution from Eq. @, and 
the N Z contribution from Eq. (35). These ef- 



fects are multiplied by the QCD correction factor 
f(a 8 ) from Eq. (22). We assume that the errors 
can be combined in quadrature. This leads to the 
net result 

S^lQ»=BGeV -> 2.62 ±0.15, B(y) =1 
2.50 ±0.17, B(y) = 2 . 

(36) 

We can compare this with the theoretical value 
Sgls — 2-62. We see that these additional terms 
improve the agreement between theory and exper- 
iment. In the approximation B(y) = 1 there is 
now excellent agreement between theory and ex- 
periment; when B(y) = 2 the experimental point 
is now below the data but still within one standard 
deviation. Since the errors are added in quadra- 
ture, the net result is a small increase in the overall 
error. 

For the purpose of completeness we will review 
the corrections that were applied to the experi- 
mental data by the CCFR group. Before solving 
for the structure functions from the neutrino and 



antineutrino cross sections (see Eq. (13) and fol- 



lowing equations), the CCFR collaboration made 
a series of corrections to the cross sections. First, 
as we have mentioned previously, the neutrino and 
antineutrino cross sections were normalized to the 
total fluxes. Next, the cross sections were multi- 
plied by four nuclear correction factors, 



norr ^iso .> ^rad w c 

<j — a x a x tr 



x a 



w 



(37) 



In Eq.(37), the term <r includes the radiative 



corrections to the cross sections, calculated from 
the prescription of Bardin and Dokochueva [4"5] , 
The term a c represents a correction for the finite 
charm quark mass since the data, particularly at 
low Q 2 , is taken in a region close to charm quark 
threshold. The term a w is a correction for the fi- 
nite W mass. The remaining correction, <t iso , was 
an attempt to account for the neutron asymme- 
try in iron. We will review this correction in some 
detail. 

On average for the CCFR target, the neutron 
excess is given by Vf = (N — Z)/A = 0.0567 
[19] . The CCFR group made an "isoscalar correc- 
tion" to the neutrino and antineutrino cross sec- 
tions. They calculated the quark and antiquark 
neutrino momentum densities on iron, via the Fi 



structure functions per nucleon for v and v on a 
non-isoscalar target, 

xq vA {x) = 2x[^(d p (x) + s p (x)) 



-A 



N 



+ -(d n (x) + s n (x))] 
— x(l — Vf)d(x) + x(l + Vf)u(x) + 2xs(x); 
xlf A (x) = x(l — Vf)u{x) + x{l + Vf)d{x); 
xq vA {x) — x(l — i/f)u(x) + x(l + i/f)d(x); 
xq vA (x) — x(l — Vf)d(x) + x(l + Vf)u{x) 
+2xs(x) . (38) 

The cross sections were then re-normalized by the 
"isoscalar correction factor," 



corr 



a(isoscalar target) 



a{v f = 0) 



cr(Fe target) a{v ! = 0.0567) 

The isoscalar correction is different for neutri- 
nos and for antineutrinos. Note that this pro- 
cess is circular - the isoscalar correction applied to 
the cross sections requires knowledge of the par- 
ton distributions, which are themselves extracted 
from the cross sections. Thus the process was ap- 
plied iteratively. An isoscalar correction was ap- 
plied to the cross sections from Eq. ( 39 ) using the 



parton distributions from Eq. (38). The structure 
functions were then determined by inserting the 
renormalized cross sections into Eq. (13). From 
the structure functions one can extract new parton 
distributions, from which a new isoscalar correc- 
tion factor could be determined. The process was 
then iterated until the difference in the extracted 
structure functions became sufficiently small. 

The isoscalar correction factor applied by the 
CCFR collaboration should account for most of 
the neutron asymmetry corrections. However, af- 
ter this correction has been applied it is then dif- 
ficult to isolate the rem aini ng contribution from 
the AF^ A term in Eq. ( fl3] ). Certainly there is a 
term present in the coupled equations, which has 
not been accounted for by the CCFR group. It 
should be straightforward to include this term in 
any re-analysis of neutrino cross section data. One 
could add this as a perturbation and could obtain 
decent estimates of the strange quark and CSV 
contributions as outlined in Sect. IIHI 

The sign and magnitude of the strange and CSV 
contributions should be similar to our estimate of 
these terms for an isoscalar target. There may 
also be small additional contributions from neu- 
tron asymmetry to xF^, which have not been ac- 
counted for by the isoscalar correction made by 
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CCFR. We have not made further corrections for 
any nuclear modification of the parton distribu- 
tions in iron. Several groups have estimated the 
magnitude of nuclear effects on parton distribution 
functions Q3HLI]. 

V. CONCLUSIONS 

The Gross-Llewellyn Smith sum rule is obtained 
from the first moment of the structure function 
CC-F3 from neutrino charged-current deep inelastic 
scattering. In principle these structure functions 
can be obtained by comparing sums and differ- 
ences of neutrino and antineutrino DIS cross sec- 
tions on an isoscalar nucleus. Previous analyses of 
the GLS sum rule have neglected potential con- 
tributions from strange quark asymmetries and 
from partonic charge symmetry violation. At the 
time, such contributions were largely unknown and 
could be assumed to be negligibly small. 

However, recently one has more quantitative re- 
sults for strange quark asymmetries from several 
groups 29 33J . All of these analyses rely on mea- 
surements of opposite-sign dimuon production in 
neutrino and antineutrino reactions on iron from 
the CCFR and NuTeV collaborations [3H|35]. All 
of these analyses obtained a positive value for 
(xs~(x)); with the exception of tnc NNPDF anal- 
ysis [3T] (which used a neural network approach, 
was relatively insensitive to the s quark distribu- 
tion and obtained very large error bars), these 
analyses found a crossover point for the s~ (x) 
PDFs that occurred at an extremely small value 
x Ri 0.01. We chose as an example the strange 
quark analysis by the NuTeV group, Mason et 
al. [55], but the correction to the GLS sum rule 
arising from the strange quark asymmetry should 
be quite similar if one used instead the CTEQ [30 , 
MSTW [3J or Alekhin [33] analyses. 

Furthermore, one now has reasonable estimates 
for contributions from valence quark CSV. First, 
there are now phenomenological analyses of parton 
distributions that include partonic CSV [5]. Sec- 
ond, there have been calculations of partonic CSV 
arising from the different electromagnetic coupling 
of photons to up and down quarks [33J 131] . Fi- 



nally, there are quark model calculations of par- 
tonic CSV [7]. We used these to estimate the 
partonic CSV contribution to the Gross-Llewellyn 
Smith sum rule. Finally, we estimated the contri- 
bution to the GLS sum rule from the fact that iron 
is a non- isoscalar target. 

The correct procedure would be to incorporate 
these effects into the initial analysis of the neu- 
trino cross sections. In particular one should take 
into account the effect of the term AF^ A (x) in 
Eq. ©. Since this term has not been included 
in previous analyses of neutrino cross sections we 
can only estimate its effect on the Gross-Llewellyn 
Smith sum rule. This has been carried out in this 
paper. To summarize our conclusions: first, the 
contributions from strange quarks, parton CSV 
and non-isoscalar effects all appear to have the 
same sign and hence to add coherently; second, 
we estimate that these effects should contribute 
an amount on the order of one to two standard 
deviations in the GLS sum rule; third, we find 
that inclusion of all of these contributions should 
bring the theoretical and experimental determina- 
tions of the GLS sum rule in agreement within 
la. In Sect. IIVDI we noted that the cross section 
corrections adopted by the CCFR group make it 
difficult to provide quantitative estimates of the 
effects of strange quarks and partonic charge sym- 
metry violation to the Gross-Llewellyn Smith sum 
rule. Nevertheless, if these corrections were to be 
integrated into a re-analysis of the neutrino cross 
sections, one should be able to obtain an accurate 
quantitative assessment of the contributions from 
these quantities. 
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